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We develop a scattering theory to investigate the multi-photon transmission in a one-dimensional 
waveguide in the presence of quantum emitters. It is based on a path integral formalism, uses 
displacement transformations, and does not require the Markov approximation. We obtain the 
full time-evolution of the global system, including the emitters and the photonic field. Our theory 
allows us to compute the transition amplitude between arbitrary initial and final states, as well 
as the S-matrix of the asymptotic in- and out- states. For the case of few incident photons in 
the waveguide, we also re-derive a generalized master equation in the Markov limit. We compare 
the predictions of the developed scattering theory and that with the Markov approximation. We 
illustrate our methods with five examples of few-photon scattering: (i) by a two-level emitter, (ii) 
in the Jaynes-Cummings model; (iii) by an array of two-level emitters; (iv) by a two-level emitter 
in the half-end waveguide; (v) by an array of atoms coupled to Rydberg levels. In the first two, 
we show the application of the scattering theory in the photon scattering by a single emitter, and 
examine the correctness of our theory with the well-known results. In the third example, we analyze 
the condition of the Markov approximation for the photon scattering in the array of emitters. In the 
forth one, we show how a quantum emitter can generate entanglement of out-going photons. Finally, 
we highlight the interplay between the phenomenon of electromagnetic-induced transparency and 
the Rydberg interaction, and show how this results in a rich variety of possibilities in the quantum 
statistics of the scattering photons. 

PACS numbers: 03.65.Nk, 42.50.-p, 11.55.-m, 72.10.Fk 


I. INTRODUCTION 

The exploration of quantum optical systems character¬ 
ized by strong photon-photon interactions has inspired a 
lot of research and extensive studies recently [ 3 - 13 . Th ose 
systems provide us with a versatile platform to investi¬ 
gate the generation and transport of non-classical light, 
as well as the behavior of single-photon sources and 
switches S Those are the basic ingredients in 

quantum-optical |l2j| and quantum information devices 

The manipulation of non-classical light typically re¬ 
quires devices displaying either strong nonlinearities Q, 
or quantum interference effects d, m- Among many 
other phenomena, they give rise to peculiar quantum sta¬ 
tistical behavior of the emitted or scattered photons, like 
anti-bunching in the generation of single photons 
or photon pairs [lH, 0. Such devices are being inves¬ 
tigated in different incarnations, including cavity QED 
[J EH, EH MM . solid state [H, ItI. f25l Eq . and circuit 
QED systems jar l30| . At the many-particle level, the 
strong interaction between the nonlinear devices and the 
photons results in the generation of many-body states, 
which can be studied in terms of dissipative versions of 
quantum spin models [3l|. In particular, the investiga¬ 
tions of the atomic steady state and the photon transmis¬ 
sion properties reveal a rich variety of quantum phases 
j32l . l33l | and photon statistics. 

In order to characterize how atomic (or any other) non¬ 
linear devices can be used to create and manipulate pho¬ 
tonic states, one can analyze the transmission spectra 
and the photon statistics, for instance, in terms of the 


second order correlation function of the photons emitted 
under the presence of weak driving light. This analysis 
is typically addressed through an input-output theory re¬ 
lating the correlation functions of the emitted photons to 
those of the atomic system in steady state E3- Those can 
be determined using a master equation approach, based 
on the Born-Markov approximation. This approach has 
proven to be very successful in most of the experimen¬ 
tally relevant situations. Even though that is a very good 
approximation for most models in quantum optics, in the 
presence of several emitters, or in certain regimes its use 
may not be justified. 

In order to analyze the transmission properties exactly, 
several elegant approaches [20|, [2l|, fUl |43I | have been de¬ 
veloped for the few photon scattering process, where the 
Born-Markov approximation is not involved. The ex¬ 
act analysis of single and two photon transmissions was 
first addressed through the Bethe ansatz approach [|20| |. 
which is equivalent to the Lippmann-Schwinger scatter¬ 
ing theory. This approach establishes the exact scatter¬ 
ing matrix (S'-matrix) between the in- and out- asymp¬ 
totic states of photons, which determines the transmis¬ 
sion spectrum and the second order correlation function 
of outgoing photons. It turns out that the Bethe ansatz 
approach is very successful in the two-photon scatter¬ 
ing by a single emitter with simple structure, e.g., the 
two-level emitter, however, the generalization to the pho¬ 
ton transmission by several emitters is difficult. The ap¬ 
proach [34] based on the input-output theory is able to 
provide the exact ^-matrix for the two photon scatter¬ 
ing by two emitters. Here, since a closed set of motion 
equations for emitter operators are required, the exact S- 
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matrix can only be obtained in some very limited cases, 
e.g., two emitters with simple structures. 

The systematic approach to the exact S'-matrix in the 
multi-photon scattering process is based on quantum 
field theory. Through either the Lehmann-Symanzik- 
Zimmerman (LSZ) reduction formalism 0 or the input- 
output theory, the exact ^S-matrix is related to the time 
ordered correlation function of emitter operators, which 
can be obtained by the functional integration of the pho¬ 
tonic bath modes. The functional integration provides an 
efficient approach to analyze the photon statistical prop¬ 
erties in the multi-photon scattering by emitters with 
complicated structures, e.g., the Jaynes-Cummings (JC) 
system [22| and the atom-coupled whispering gallery res¬ 
onator 24J. However, the exact ^-matrix is only able 
to describe the evolution of photonic bath in the asymp¬ 
totic limit, thus, it fails to depict the transient dynamics 
of emitters, e.g., the single-photon detection in the super¬ 
conducting edge sensor (4J-|46|. Recently, a generalized 
input-output theory [36| was established by the subtle 
combination of the input-output theory and the quan¬ 
tum regression theorem, which allowed us to investigate 
both the transient dynamics and the photon statistics in 
the asymptotic limit. Here, a generalized master equa¬ 
tion is obtained to describe the transient dynamics of 
several emitters under the presence of few incident pho¬ 
tons. The generalized input-output theory and master 
equation involve the Markov approximation, thus, they 
are not able to provide the exact results. The use of the 
Markov approximation in the presence of several emitters 
need to be justified. 

In this paper, based on the path integral formalism, 
we develop the scattering theory to characterize the full 
time evolution of the global system (including the emit¬ 
ters and the photonic bath) exactly for the few-photon 
scattering by several emitters with complicated struc¬ 
tures. Our theory provides the transition amplitudes 
from the arbitrary initial state to the corresponding final 
state without the Markov approximation, where both the 
quantum statistics of scattering photons and the tran¬ 
sient dynamics of emitters can be analyzed exactly. In 
the Markovian limit, the exact results from our theory 
perfectly agree with those from the quantum regression 
approach, which justifies the validity of the Markov ap¬ 
proximation. Similarly, the generalized master equation 
is re-produced to describe the dynamics of emitters in 
the presence of few incident photons. Here, we empha¬ 
size that the generalized master equation establishes the 
close relation between the properties of photons emitted 
under the present of weak driving light and those in the 
few-photon scattering by the emitters. 

Using the two paradigmatic examples, i.e., the few- 
photon transmission to the two-level emitter and the JC 
system, we examine the correctness of our theory in the 
photon scattering by the single emitter. Here, the de¬ 
veloped scattering theory re-produces the well-known re¬ 
sults PqI-ElI [47) for the transmission spectra and the sec¬ 
ond order correlation function. For the few-photon scat¬ 


tering by several emitters, we investigate the condition 
of the Markov approximation in detail. In particular, we 
show that the Markov approximation is valid under cer¬ 
tain conditions relating the bandwidth of the dynamics 
and the distance of separation between emitters. In the 
non-Markovian regime, the exact results exhibit peculiar 
features associated with retardation of pulses between 
emitters. 

By the developed scattering theory, we explore some 
novel phenomena in two new situations. For the photon 
scattering by a single two-level emitter in the half-end 
waveguide, we show how the two-level emitter can gener¬ 
ate entanglement of out-going photons. For the photon 
transmission in the array of atoms under conditions of 
electromagnetically induced transparency (EIT) [48[ and 
coupled to Rydberg levels, we highlight the EIT phe¬ 
nomenon and the Rydberg interaction. Here, the second 
order correlation functions of emitted photons, the co¬ 
propagation of dark polaritons and the collision of the 
counter-propagating polaritons in the transient process 
are analyzed, which show a rich variety of the quantum 
statistics of photons and Rydberg excitations. The devel¬ 
oped scattering theory is proven to be a very efficient and 
systematical approach to investigate the quantum statis¬ 
tics of photons in the array of the interacting emitters 
with complicated structures. 

The paper is organized as follows. In Sec. II, the five 
models are introduced, which are the two-level emitter 
and the JC system coupled to the photonic waveguide, 
the two-level emitter in the half-end waveguide, and the 
photon scattering by an array of two-level atoms and EIT 
atoms coupled to Rydberg-levels. In Sec. Ill, the de¬ 
veloped scattering theory is established, where the exact 
S'-matrix relates the asymptotic incident and final states. 
For the transient process, the generalized master equa¬ 
tion is derived to characterize the time evolution of the 
quantum emitters. In Sec. IV, two paradigmatic exam¬ 
ples, i.e., the photon transmission to the two-level emitter 
and the JC system, are used to examine our theory and 
illustrate how our method works in the photon scattering 
by the single emitter. In Sec. V, we analyze the validity 
of the Markov approximation for the photon scattering 
by several emitters. In Sec. VI, we show that the gen¬ 
eration of entanglement between two scattering photons 
by a single two-level emitter in the half-end waveguide. 
In Sec. VII, the photon transmission to an array of EIT 
atoms coupled to Rdyberg levels is investigated. In Sec. 
VIII, the results are summarized with the outlook. 


II. SCATTERING MODELS 

In this section, we introduce the models we are going to 
use in order to investigate the transmission of waveguide 
photons interacting with quantum systems. Those could 
be sing le emitters (e.g., a multi-level atom jji fill. I20I. EL 
Eol Isd ] . an atom coupled to a cavity mode [it, i;2| ), or an 
array of emitters coupled, for instance, to Rydberg levels 
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The model Hamiltonian has the form H = H w + H sys + 
Hr + Hr f and contains four parts: (a) The free propa¬ 
gation of photons in the waveguide is described by the 
Hamiltonian 

= y^[(& 0 + k)r\r k + (fco - fc)Zjjfc], (1) 

k 

where (r£) and Ik (l\.) are the annihilation (creation) 
operators for the right- and left- moving modes with the 
central frequency fco in the waveguide; (b) The system 
Hamiltonian H sys describes the quantum emitters, and 
will be explicitly given latter for different examples; (c) 
The interaction between the waveguide photons and the 
emitters is given by the Hamiltonian 

Hr = J2lVHi r H x i) + Vr^l\xi)}Oi + H.c., (2) 
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where the left- and right-moving photon fields r(xj) = 
E k r k e i{ - k+k ^ Xi /VL and l{ Xi ) = J2k l ke i{k ~ ko)xi /VL 
couple to the operator Oi of the i-th emitter at the po¬ 
sition Xi with coupling strengths T^ i . In momentum 
space, the interaction Hamiltonian reads 

Hr = J2( r k°^+ + llOk,-) + ^c., (3) 

k 

where the collective operators Ok,± = 
Y)i VT ~iOie ~ l ^ k±k ^> Xi and we focus on the symmet¬ 
ric coupling case, i.e., T r ^ = T^ = T^. Without loss of 
generality, we choose the position of the first emitter 
i m. 1 at the origin x\ =0; (d) Apart from the decay to 
the waveguide, there may exist other decay channels to 
free space (spontaneous emission outside the waveguide). 
We describe this through the Hamiltonian Hr f . Here, 
each emitter couples to some system operator Oi with a 
JC-type coupling ©• 

We now present some relevant examples for the system 
Hamiltonian H sys , and which will be analyzed in detail in 
the following sections. As shown in Fig. QJi, the simplest 
example consists of a two-level system with the energy 
level spacing uj e = fco between states \e) and | g). The 
Hamiltonian is H sys = uj e \e) (e| and O l= \ = O, = O = 
o \g)(e\. 

The second example is the JC system shown in Fig. 
nr The Hamiltonian 

H sys = w c a f a + u e |e) (e| + g(a f | g) (e\ + H.c.) (4) 

describes a single two-level system, with the energy level 
spacing cj e , coupled to a cavity mode of frequency uj c . 
The annihilation operator a corresponds to the cavity 
mode, and the coupling constant is g. Here, the cavity 
mode directly couples to the waveguide, i.e., 0*=i = O = 
a. The excitations of cavity field and the two-level system 
can decay into free space; thus, Oi would be a and | g) (e\. 


FIG. 1: (Color Online) Four examples of the problem studied 
here: (a) Single two-level system coupled to a waveguide; (b) 
JC system coupled to a waveguide; (c) Single two-level sys¬ 
tem coupled to a one-sided waveguide; (d) An array of atoms 
coupled to a waveguide. 


We analyze the validity of the Markov approximation 
and the retardation effect using the third example. Here, 
as shown in Fig. [Hi, the waveguide photons couple to an 
array of two-level emitters with frequency uj e = fco and 
lattice spacing d. The emitter Hamiltonian is 

H sys = + ±U 0 b\btbibi), (5) 


where the hardcore boson bi is introduced to describe a 
two-level emitter, and the hardcore behavior is obtained 
in the limit Uo oo. The wave guide photons couple 
to the collective modes Ok,± = %/F/L Yi e~ l ^ k±k °^ Xi bi 1 
where the decay rates T are taken to be the constant 
r= T. The excitation decays into the free space, 
i.e., Oi = bi. 

The fourth example is a single emitter in front of a 
mirror, as shown in Fig. []J. The system Hamiltonian is 
H sys = ^emitter + and the collective operators are 


Ok,+ 


Ok- 


Q e -i(k+k 0 )x 0 

Q e -i(k-k 0 )x 0 



(6) 


where the boson mode b is introduced to describe the 
mirror in the limit T^ —>• oo. The decay of the emitter 
excitation into free space is characterized by Oi = \g) (e|. 

The fifth example is an array of EIT atoms coupled to 
Rydberg levels, with lattice spacing d, as it is schemati- 
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cally shown in Fig. [TJi . The emitter Hamiltonian is 

H s ys = ^[cj e et ei + cjgstsi + + H.c.)] 

+JT H c + 2 ( 7 ) 

ij 

Here, we use hardcore bosons e and s to describe a 
Rydberg-EIT atom. There are three atomic levels: the 
ground state | g), the excited state \e) with frequency u; e , 
and the Rydberg state | s) with frequency uj s . The hard¬ 
core behavior is obtained through 

#HC = -y + s i s i)( e i e i + 4 s * - !) (8) 

i 

in the limit Uo oo, which projects out the states 
|ee), |es), and \ss) of each atom with double occupa¬ 
tions. The atoms in the Rydberg state | s) experience 
long-range interactions, Uij, (52]. For instance, we can 
take the Van-der-Waals interaction Uij = Cq/ \xi — xj | 6 , 
the dipolar interaction Uij = C 3 J \xi — Xj\ , or the uni¬ 
form interaction Uij = Co [36|. The transition be¬ 
tween states \e) and |s) is induced by a classical field 
of frequency uJd and corresponding Rabi frequency fi. 
The waveguide photons couple to the collective modes 
Ok,± = v^T JT e~ l ^ k±ko ^ Xi ei/y/~L, where the decay rates 
r r p)^ are taken to be the constant T r (i),i = T. The e- 
excitations can decay into the free space, i.e., Oi = e*. 

For convenience, we transform the Hamiltonian in the 
rotating frame rk —» rke~ lkot and Ik lke~ lkot . The sys¬ 
tem operators are also transformed correspondingly such 
that the total Hamiltonian is time-independent. As a re¬ 
sult, in the rotating frame, H w = J2k k(r\rk — /£/&), and 
the frequencies in the system Hamiltonian are replaced 
by the detunings. This will be used in the five examples 
presented in the following sections. 


III. GENERAL FORMALISM 

In this section, we show the general formalism to study 
the full dynamics of the global system during the few pho¬ 
ton scattering process. In Secs. IIIA and IIIB, we briefly 
review the results based on the Markov approximation in 
Ref. 136]. Here, a different method, i.e., the displacement 
transformation, is introduced to relate the transition am¬ 
plitude between arbitrary initial and finial states with the 
correlators of emitter operators. Based on the quantum 
regression theorem, these correlators are obtained by the 
effective Hamiltonian of the emitters, where the quantum 
regression theorem requires the Markov approximation. 
In order to obtain the exact result and examine the valid¬ 
ity of the Markov approximation, in Sec. IIIC we use the 
path integral formalism to derive the exact transition am¬ 
plitude, where the non-Markovian effects are taken into 
account. 


In both approaches, the S'-matrix of the asymptotic in- 
and out- states and the transition amplitudes in the tran¬ 
sient process are obtained, which agree with each other in 
the Markov limit. It turns out that for the single emitter 
coupled to the photon with linear dispersion, the Marko¬ 
vian result is proven to be exact by the path integral 
approach. For the photon scattering by several emitters, 
the validity of the Markov approximation is analyzed in 
Sec. V in detail. In the Markov limit, both approaches 
give rise to the generalized master equation |36j governing 
the dynamics of the emitters during the scattering pro¬ 
cess, where the generalized master equation establishes a 
close relation between the transient behaviors of emitters 
in the presence of few incident photons and the dynamics 
of emitters under the weak driving light. 


A. S- matrix by quantum regression theorem 

In this section, based on the quantum regression the¬ 
orem, the displacement transformation is used to derive 
the transition amplitude between arbitrary initial and 
final states, where the quantum regression theorem in¬ 
volves the Markov approximation. 

We introduce the transition amplitude 

A{T) = sys (v? out | (out| (*/-*<) |i n ) |^ in ) sys (9) 

from the initial state |in) |<A„) sys = |0) free |V’in) w |<An) sys 
to the final state |out) |(/?out) S y S = 10) free \^P OUt) w |^OUt) S y S 
during the time T = tf—ti. At the instant the waveg¬ 
uide photons are in the state |t/>i n ) w , and the initial state 

of the emitters is |^ in ) sys = 7iL |G) 8ys , where |G) sys and 

y/ n denote the ground state and some creation operator 
of the emitters, respectively. Similarly, at the instant tf, 
the corresponding final states are |^ 0 ut) sys = Tout |G) sys 
and |V ; out) w? respectively. Here, we focus on the transi¬ 
tion process without excitations leaking to the free space, 
and the initial and final states of the free space are the 
vacuum state |0) free . 

The initial and final states of waveguide photons can 
be generally written as 

IV>in) w = E V’in({nfca})RKa), 

{nkcx} k 

IV ; out) w = ^ ^ ^t(K a })nK), (10) 

{m ka } k 

where {n ka } = {n klCt ,n k2a ,...} and {m ka } = 
(fUfcicn TUfecn •••} are the number distribution of photons 
with different momenta ki in the initial and final states, 
respectively, and a = r, / denote the right- and left- 
moving modes. We notice that the relation 

1 Q n kcx 

\n ka ) « lim ———^\J ka ) (11) 

J ka ^o y/n ka \ oJ k a 
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between the Fock state and the unnormalized coherent 
state I Jka) = Y,n k J ka \ n ka) /leads to the coher¬ 
ent representation of the initial and final states 


where 

‘Fin 

F’out 


|V ; in) w — Fin | {Jka}) 5 
Iout) w = F3 U t | {Jka}) 5 ( 12 ) 


lim 

{Jka}—> 0 


{rika} ka 


1 () nfcQ; 


(13) 


lim n E V’out({m fca })TT 

{•W->o ^ -O 


1 ^TYlka. 

y/mkJ 5J^ a ' 


In terms of Eq. m , the transition amplitude reads 


A(T) = (14) 


It follows from Eq. m that the functional deriva¬ 
tive of Aj(T) determines the transition amplitude A(T), 
which is composed of Fourier transforms of time ordered 
correlators 

G(T) = (l~ 7out (tf)Oi(ti)...O n (t n )7hi(^*)) (20) 

on the ground state |{0^}) b |G) . Here, Oj(tj ) = 

W (t)OjU(t) is given by the emitter operators Oj = 0/c J5 ± 

and °l,±- 

By the quantum regression theorem, as shown in Ref. 
[33, the bath degree of freedom can be traced out and 
the correlator 

G(T) = sys (G\ F7out(^/)O e ff,l(H)---O e ff,n(^n)7 i n(^) l^) S ys 


where the transition amplitude 

Aj(T) = sys (^outlb {{Jka }| 6 |{Fa})b | ( ^in) sys , 

(15) 

and |{J/ca})b = 10)^ | {Jka})- The transition amplitude 
m can be evaluated by either the quantum regression 
theorem or the path integral approach. 

We first show the result from the quantum regres¬ 
sion theorem, where the Markov approximation is re¬ 
quired. In the interacting picture, the time evolution 
operator hi = T exp [—i f* f dtH(t)\ is determined by the 
Hamiltonian H(t) = H sys + Hr f + Hr(t), where T 
is the time-ordering operator and the interaction part 
H v {t) = is 

H r (t) = E(^O fe ,+e ife4 + l{O k ^e~ ikt ) + H.c.. (16) 

k 

In terms of hi , the amplitude (fl5l) reads 
Aj(T) = sys ((/?out|b ({Jka,out}\M \{ Jka,in}) |^in) sys j 

(17) 

where Fa,in = Jka^ l£kcxti , Jka, out = Jkae l£k ° itf , the dis¬ 
persion relations £k,a = &ak, and a r p) = d= for the right- 
and left- moving modes, respectively. 

The displacement transformation U : | {Jka}) = 

e Efcal J fc«l / 2 U \{0ka}) is introduced to rewrite the am¬ 
plitude 

Aj{T) = exp( E \JkA 2 e ~ i£k ’ aT ) x (18) 

k,a=r,l 

sys {^P out |b ({0/ca}| U |{0/ca})b I ( ^ > in) sys 

where the state of the waveguide photon is transformed 
to the vacuum state, and the time evolution operator 

U = dt i H ^A H d Ml determined by H{t) and 

the driving term 

H d (t) = E (19) 

k,a=r,l 

+Jk, a Ol aa e- i ^ t ~ ti \ 


becomes the average value of emitter operators 
Oes,j(t) = Ul s {t)OjU e s{t) on the ground state |G) syg . 
The time-evolution operator U e s(t) = exp(— iH e ^t) is 
given by the non-Hermitian effective Hamiltonian H e ^ = 
H sys + i^decay, where 

H ,decay = f^Oi *7 

i ij 

( 21 ) 

and T f : i is decay rate to the free space of the emitter at 
the position X{. 

Based on the quantum regression theorem, each term 
in A(T) is related to the correlator of emitter opera¬ 
tors O e fij(t) governed by the effective Hamiltonian H e ff. 
The re-summation of these terms results in A(T) = 
F* ut Fin^4j(F) with the compact form 

Aj(T) = ByB (¥»out|Z4ff|v>in) ByB e^.<*=r. 

( 22 ) 

where U e a = Te~ 1 ^ dt ( H e«+H d ) ^ p] le equations (TT~TT) 
and ([22] ) establish the relation between the transition am¬ 
plitude A(T) and the time ordered correlators of emitters 
governed by the effective Hamiltonian H e ^. 

By different choices of 7 in , 7 out , U, £/, F in and F ou t, 
the dynamics of the global system can be fully character¬ 
ized. Hereafter, we refer the choice of 7 in(out ) 5 ij )> and 
Fi n (out) as the boundary condition. For instance, A(T) 
with the boundary condition 7 i n = 7 out = I (identity 
operator), U —> — oo, and tf —> oc leads to the photonic 
F-matrix. For the boundary condition 7 i n = /, 7 out 7 ^ /, 
U = 0, and tf = T, A(T) describes how the incident 
photons transform to the emitter excitations and propa¬ 
gate in the transient regime. For the boundary condition 
7in 7 ^ 1 , 7out = /, ti = 0, and tf m T, the behaviors 
of spontaneous and stimulated emissions can be investi¬ 
gated by A(T). Using Eqs. (fT4]h (|2T]h and ([22]) . together 
with different boundary conditions, we shall study the 
few photon scattering process in the five models in Secs. 
IV-VII. 
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B. Generalized master equation 

In order to study the dynamics of emitters during the 
scattering, one can either use the transition amplitude 
(HI with proper boundary conditions or derive the mas¬ 
ter equation of emitter reduced density operator by trac¬ 
ing out the photonic degree of freedom. In quantum op¬ 
tics, the initial states of the photonic bath are usually 
the vacuum state, the thermal state, and Gaussian states, 
where the conventional master equation is obtained based 
on the Born-Markov approximation. 

During the scattering process, the initial state of pho¬ 
tons is dramatically changed. For instance, the single 
incident photon resonant with the transition energy of 
the two-level emitter is totally reflected, as predicted by 
the single photon scattering theory [3|, Uj], [2CJ, [2lj , where 
the photon in the initial asymptotic state r k |0) is to¬ 
tally scattered to the photon in the out-going asymptotic 
state l^_ k |0). However, the conventional master equation 
assumes that the initial state of the photonic bath is un¬ 
changed, and the state of the global system is always the 
product state of the emitter and the photonic bath E3- 
This assumption contradicts with the exact result of the 
single photon scattering theory, thus, the conventional 
master equation breaks down. The generalized master 
equation is required to investigate the dynamics of emit¬ 
ters in the few-photon scattering process. 

For the initial state of few photons in the waveguide, 
the evolution of the system state is described by the re¬ 
duced density matrix p s (T) = Trb a th[Wp(0)Wl], where in 
terms of the coherent state \{Jka}) 

P( 0) = ^in[Psys(0) 0 \{Jka}) h ({Jka} |]- (23) 

The displacement transformation U relates the density 
matrix 

Ps(T) = T* n T in eZ^ J ^ 2 pj (T) (24) 

to the generating density matrix 

PJ (T) = Tr hath [Ujp sys ( 0) ® |{0 feQ }> ({0 ka }\U% (25) 

where Uj = U^UU describes the evolution of emitters 
under the driving field. In pj(T), the initial state of 
photons is transformed to the vacuum state, thus, the 
conventional master equation can be used to describe 
the time evolution of pj(T). Under the Markov approx¬ 
imation, the master equation governing the evolution of 
Pj(T) reads 

d T pj(T) = -i[H s ys , pj(T )] + Cpj(T ), (26) 

where the initial condition is pj(0) = p sys (0). The dis¬ 
placement transformation induces the driving term in 

H S ys — H S y S T ^ ^ y/YiYjOi Oj sin(&o |%i %j |) 

*3 

+ y E( J l*O k ,a a e ia “ kT + Jk, a Ot taa e-^ kT ), (27) 

k,a 


and the Lindblad operator is 

Cpj(T) = 2 y] Oipj(T)0] cos[fc 0 (a;j - xj)\ 

ij 

- J 2 v/nFcos[*d(*i - (28) 

ij 

+ Y J TfAzd i pj{T)d\ - {ojOi, P j(T)}]. 


The generalized master Eqs. ([24]) and ([26]) are the 
main results in this section, which lead to some signif¬ 
icant results: (a) In principle, once pj(T ) is obtained, 
p s (T) is determined through Eq. ([24]) , which describes 
the transient dynamics of emitters during the scattering 
processes. In practice, the evolution of emitters for the 
few incident photons can be studied by the perturbative 
expansion of the classical sources Jj . C?Q , and J k a . For in¬ 
stance, for the single incident photon, Eq. (1241) contains 
at most the second order derivative S 2 pj(T)/SJ* a ,SJk a , 
the second order perturbative expansion of Jk, a and a 
in pj(T) determines p s {T). Here, the zero order and 
second order contributions lead to the emitter reduced 
density matrix for the single incident photon. The zero 
order contribution p sy s(0) describes the unchanged emit¬ 
ter state without interacting with the photon, while the 
second order contribution Pj(T) describes the dy¬ 

namics of emitters responding to the single incident pho¬ 
ton. Similarly, for two incident photons, the reduced den¬ 
sity matrix can be determined by the forth order expan¬ 
sion of J k ,a and J k a . 

(b) In quantum optics, we are interested on the trans¬ 

mission properties of weak probe light to emitters, i.e., 
the transmission spectrum and the quantum statistics. 
We notice that Eq. (l26j) just describes the emitters un¬ 
der the classical probe light with the strength Jk a • If the 
driving field is weak, we can solve Eq. ([26]) by the per¬ 
turbative expansion of Jk, a and As we discussed 

above, the second and forth order perturbative expan¬ 
sions describe the transient dynamics of emitters in the 
single and two photon scattering processes. As a result, 
the quantum properties of out-going photons and emit¬ 
ters under the weak driving light can be studied by the 
few photon scattering theory, where only the time evo¬ 
lution in the forward path is involved, as shown in Eq. 
(|22|) . We shall show in Sec. VII that for the emitters 
with complicated structures the scattering theory is still 
able to provide the analytic results, which perfectly agree 
with the numerical results from the master equation ap¬ 
proach. The few photon scattering theory enables us to 
analytically investigate the quantum properties of pho¬ 
tons from emitters under the weak diving field. 

(c) If we set the external source to be zero, Eq. ([26]) 
agrees with the master equation [3.1] for the photonic 
bath initially in the ground state. 
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C. Exact ^-matrix by path integral approach 

The equations (fT4]h (|2T|) , and ([22] ) in Sec. IIIA relate 
the transition amplitude to the correlators of emitters 
governing by the effective Hamiltonian, where under the 
Markov approximation the quantum regression theorem 
leads to the instantaneous effective Hamiltonian. In this 
section, we show the exact result of the transition ampli¬ 
tude d, from which the result obtained by the quantum 
regression theorem is proven to exact for the single emit¬ 
ter coupled to the waveguide. The validity of the Markov 
approximation and some non-Markov effects for several 
emitters coupled to the waveguide will be investigated by 
the exact transition amplitude in Sec. VI. 

Following the procedure of path integral formalism, 
we discretize the evolution time T by N oo instants 
{£i....£jv} and insert coherent state basis at each instant 
U. By integrating out the free space modes we obtain 

Aj{T) = I Dlsystembo^/b^)^ 

x J D[a k ,a* k ]e^ J ^ ak{t ^e iS , (29) 

where the action 

S = I dt{Y\r* k {idt-k)r k +r k {id t +k)l k ]-H T } (30) 

Jt * k 

describes free propagation of the right- and left- moving 
photons in the waveguide and the interaction with the 
emitters, and f D [system] is the integral over the emitter 
field. The action of the emitter is 

s sys = s£l + i f tf dt V r f ,id]di, (31) 

Ju i 

where the first term is the action of emitters, and the 
second term describes the decay to the free space. 

By SS/Srl = 0 and SS/Sl^ = 0, the classical motion 
equations read 

(idt ~ k)r k> ci - O k ,+ = 0, 

(idt + k)l k , d - O k) - = 0, (32) 

which give the classical paths 

r k ,c\(t) = Jk,re- lk(t ~ u) - i f dsO k!+ (s)e~ lk( ' t ~ s \ 

Jt, 

h,At) = Jk,ie iHt - u) - i J dsO k ,-(s)e ik ^- s \ (33) 

Following the saddle point method, we expand the pho¬ 
ton fields rk = He,cl + Srk and Ik = h,ci + Slk by the 
quantum fluctuation fields Srk and 81 k around the clas¬ 
sical paths, and integrate out the fluctuation fields in 
Aj(T). Due to the quadratic form of the action ([30]) . 

Aj(T) = exp(^|J fe , a | 2 e-^ T ) (34) 

k,a 

X J D[system}'j out (tf)j* n (t i )e lSet{+lSj 


is obtained exactly, where the nonlocal effective action 
S e ff = S sys + S re of emitters is determined by 

S ve = i f dt f ds^2\Ol t+ {t)O k ,+{s)e~ lk{t ~ s) 

Jt, Ju k 

+Ol_(t)O k ,.(s)e ik ^-% (35) 

and the source terms Sj = — f* f dtH&(t). 

The equations (IM1) and ([35] ) are the central results in 
this section. It follows from Eq. m that the functional 
derivatives of Aj(T) lead to the transition amplitude 
A(T). We also notice that, apart from the phase factor, 
Aj(T ) is the generating functional [53] of the emitter cor¬ 
relators, namely, the functional derivatives of Aj(T ) give 
the correlation functions of emitter operators. As a re¬ 
sult, the transition amplitude A(T) is determined by the 
correlation functions of emitter fields, where the effective 
action S e s is generally time nonlocal. 

For the single emitter coupled to the waveguide with 
the linear dispersion, the action 5 e ff turns out to be 
time-local, and the effective Hamiltonian H e q is deduced. 
Here, the renormalized part has the time-local form 

S re = iT J f dtO*(t)0(t), (36) 

and effective Hamiltonian 

iJeff = H sys -iTfd ] d- iTO ] O (37) 

of emitters is obtained. In fact, the functional integral in 
Eq. (O is just the amplitude sys (v? Q ut| Weff l¥ , in) sys in Eq. 
(j55J), thus, the result based on the quantum regression 
theorem agrees with the exact result (1531) . 

For several emitters, the renormalization part 

S re = i f f dt'S'^W j O*(t)O j (t-\x i - X j\)e ik °^ Xi ~ x ^. 

Jt < a 

is time nonlocal. If the wavelength of frequency fluc¬ 
tuations around fco is much larger than the size Nd of 
the emitter array, the time delay effect can be neglected, 
which leads to the time-local action 

Sre = i f f dtY J VTiT~jO*(t)O j {t)e iko|a '*“^ 1 , (38) 

Jt > a 

and the effective Hamiltonian (I2T1) . In this limit, the 
result ([22]) agrees with the exact result (IM1) for several 
emitters. 

In the end of this section, we briefly summarize the 
results. We use two approaches to obtain the transition 
amplitude and the generalized master equation for the 
photon scattering process, where the derivation of the 
generalized master equation ([26]) from the path integral 
approach is left in the Appendix. The relation between 
the few photon scattering by emitters and the emitters 
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under weak driving light is established. For the single 
emitter case, the result ([22]) based on the quantum re¬ 
gression theorem is proven to be exact. For several emit¬ 
ters, in the Markovian limit, the result (|22|l agrees with 
the exact result from the path integral approach. 


IV. PARADIGMATIC EXAMPLES 

In this section, as paradigmatic examples, the few- 
photon scattering by the single emitter is studied by 
the scattering theory developed in Sec. III. The emit¬ 
ter is chosen to be the two-level emitter or the JC sys¬ 
tem, where the decay to the free space is neglected, i.e., 
Tf = 0. For the two-level emitter, we investigate the 
transmission spectra and the quantum statistics of out¬ 
going photons for the single and two incident photons. 
For the emitter initially in the excited state without the 
incident photon, the spontaneous emission of the emit¬ 
ter is studied. The transient process, i.e., the response 
of the two-level emitter in the ground state to the single 
photon wave-packet, is analyzed by the exact transition 
amplitude of the emitter in the excited state. For the 
stimulated emission of the two-level emitter initially in 
the excited state, the wave-packet shape of the single in¬ 
cident photon is designed, such that the probability of 
stimulated emission is maximal. For the JC system, we 
show the transmission spectra and the quantum statistics 
of scattering photons. By these two examples, we exam¬ 
ine our theory by_the comparison with the well-known 
results in Refs. [20M22l. l47j | . which justifies the correct¬ 
ness of the developed scattering theory. 


A. Two-level emitter 

For the two-level emitter, O = a~ and the effective 
Hamiltonian H e ^ = H sys — follows from Eq. 

{37D, where H sys vanishes in the rotating frame since the 
transition frequency uo e = ko is resonant with the central 
frequency of the waveguide. It follows from Eq. (fl9l) that 
the external source term is 

Sj = ~\j\ [ h dt 

^ k,a=r,l 

+Jk,*o + {t)e- iaak{t - u \ (39) 

For the scattering of single right-moving photon with 
the momentum k by the emitter in the ground state, the 
boundary condition in the asymptotic limit U —» — oo is 
7 i n = I and T \ n = lmxrj fc }_^ 0 d/dJk,r- By the boundary 
condition 7 ou t = I and J'out = lim(j }_,. 0 6/6J Pt ( r j) in 
the asymptotic limit tf —> oo, Eqs. (THl) and (pHl) lead to 
the reflection and transmission coefficients 

-i r 

Rk&p, — k 1 ' 7^dp, — ki (40) 

/CO + H 


and 


T k S Ptk = (41) 

in the asymptotic limits U — oo and tf oo. For the 
scattering of two right-moving photons with momenta k\ 
and k 2 by the emitter in the ground state, the boundary 
condition in the asymptotic limit — oo is 7 i n = I and 
T \ n = lim.rj fc }_^Q S 2 /5Jk 1: rfiJk 2 ,r- By the boundary con- 
dition 7out = / and Jy.ut = S 2 /SJ- Pu i6J- P2 f 

in the asymptotic limit tf — oo, it follows from Eqs. (fl4l) 
and (IM1) that the ^-matrix 

q 7 7 _ p 2 / + °° dt' 1 dt' 2 dtidt 2 ipit ^ +ip2t ' o _ ikltl _ ik2t2 

* PlP2;fclfe2 “ J_ x (2 tt ) 2 

(Ta~(t' 2 )a + (h)a + (t 2 )) , (42) 

of two reflected photons with momenta —p\ and — p 2 is 
given by the time ordered correlator of spin operators. 
Here, the time evolution cr^it) = a ± e ±rt governed by 
iJ e ff determines the four-point correlator, which leads to 
the S'-matrix 


&piP2',kik2 — Rk\Rk 2 (dp\k\dp2k2 d" dp 1 k 2 dp 2 k\) 
iV 2 ^p 1 +p 2 ,fe 1 +fe 2 l k 2 -\- 2iT) /,q\ 

tt (pi + ir)(p 2 + ir)(A; 1 + ir)(A: 2 + ir)' 1 J 

The wavefunction of two reflected photons is obtained 
directly from the Fourier transform of Eq. (1431) [U] as 

tp(x 1 ,X 2 ) = \ J d Pi d P 2 e WlXl+W2X2 S PlP2 . klk2 

= e iEx ‘R kl R k2 [cos(kx ) - e& B - r M](U) 

where the total momentum E = k\ -f k 2 , the relative 
momentum k = {k\ — k 2 )/2, the center of mass coordinate 
x c = (xi+x 2 )/2, and the relative momentum x = x\ — x 2 . 
The results (1451) and (144j) agree with those in Ref. j|20|, Hl| . 

For two photons with the resonant frequency k\ = k 2 = 
0, the reflection coefficients are R/ Cl = Rk 2 = 1 and the 
wavefunction 


^(r c , x) = 1 — e r|x| . (45) 

In Sec. Ill, we prove that the quantum statistics of out¬ 
going photons from the emitter under the weak driv¬ 
ing light can be characterized by the two photon wave- 
function xp(x c ,x). More precisely, for the weak driving 
light with the zero detuning — cj e , the sec¬ 

ond order correlation function g (y2 \x) = \^{x c ^x)\ 2 of 
the emitted photons can be obtained by the two photon 
wavefunction vp(x Cl x) from the scattering theory. Here, 
g^ 2 \x) = (1 — e -r ^l) 2 display the anti-bunching behav¬ 
ior. 

The spontaneous decay of the two-level emitter in the 
excited state can be investigated by the boundary con¬ 
ditions 7 i n = cr _ , Jin = 1 and 7 ou t = /, ^out = 
limpq^o at the initial and the final instants 
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X A o ' 

(a) (b) 


FIG. 2 : (Color Online) Emission and absorption of a single 
photon by a single emitter: (a) The spatial amplitude |'0(x )| 2 
of a right-propagating single photon wavepacket, produced 
by an emitter initially prepared in the excited state at T — 
0. The wavepacket shown here is after an evolution time 
T = 10/r. (b) The probability \A(T )\ 2 of the emitter in the 
excited state for the incident wavepacket with the width I /7 
and localized at xq. Here, the emission rate T = 1 into the 
waveguide is taken to be the unit. 


U = 0 and tf =T. The equations (fl4l) and ([34]) give the 
amplitude 


A(T) = 



—ia a pT _ g —FT 

cr a p + iT 


(46) 


to detect single out-going photon with the momentum p. 
The Fourier transform gives the wavefunction 


^(*) = J ALe^A(T) 

= —iy/fe~ r ^ T ~ aaX ^9(T - a a x)0(a a x) (47) 


in the coordinate space for the right- and left- mov¬ 
ing photon. Since the emitter couples to the left and 
right moving modes symmetrically, the wavepackets of 
the right and left moving photons are symmetric with 
respect to the origin. As shown in Fig. [ 2 ji, \^p(x)\ 2 for 
the right moving photon exhibits the Lorentzian shape 
with the different widths 1/T. 

The reverse process of the spontaneous decay is 
the response of the two-level emitter in the ground 
state to the single photon wavepacket f{ n (k). Here, 
the boundary conditions are 7 i n = /, T\ n = 
lim.{j fc }— >.0 J dkfi n (k)5/SJ^^ r and 7out = & and ;>ut = 1 
at the instants t% — 0 and tf = T, respectively. The 
equations ([14]) and (l34|) result in the amplitude 


A{T) = J dkf in (k) J dte lkt (a ( T)cr + (t )) 

= ^ J dkf in (k) 


e~ ikT - e~ TT 


k + iT 


of the emitter in the excited state. For the right-moving 
photon wavepacket 


7 e 


-ikx 0 


7T k + 27 


with the width I /7 initially at the position xo < 0 , the 
residue theorem gives 


A(T) = ^A[ e -n^o+T) _ e —7(*o +T)]0( Xo + T ). (50) 
r - 7 

In Fig. [2Jd, we show |^4(T )| 2 of the emitter in the excited 
state for different widths 7 . 

In the stimulated emission, the two-level emitter is ini¬ 
tially prepared in the excited state at the instant G = 0 , 
and the wavepacket of single right-moving photon is de¬ 
signed to realize the maximal probability of emitting two 
right-moving photons in the asymptotic limit tf 00 . 

With the initial and final boundary conditions 7 i n = 
c“, Jin = limjjjj-Mj / dkf in (k)6/6J kt r and 7 out = I 
Tout = lim{j fe }^o d 2 /5J plt rfiJp2,r, Eqs. (JH]) and (JM]) re- 
suit in the ^S-matrix 


S, 


'piP2,k ~ 


+ 


h 

r 


J dkf in (k)[S kp 


P2 + iT 


1 


1 


1 


2tt pi + iT pi — k + i0+ p\ + P 2 — k + zT 
+(pi^p 2 ) (51) 


of two right-moving photons with momenta pi and P 2 . 
The Fourier transform of S PlP2i k results in the wavefunc¬ 
tion 


^(x u x 2 ) = J ^l SpW 2 , k e i ^ +i ^ (52) 
= J dxf in (x)[B(x 1 ,x 2 ;x) + B(x 2 ,x k ;x)], 


in the coordinate space, where 


B[x\,X 2 \x) = —iVTe rx2 0(—x 2 )[S(x — xi) 

-Te~ T{x - Xl) e{x 2 - x)6(x - ®i)],(53) 


and 


/ rlk 

—Uk)e ikx (54) 

describes the incident wavepacket in the coordinate 
space. 

We design the shape of the wavepacket fi n {x) to max¬ 
imize the stimulate emission probability 

Pst = ^ J dxidx 2 \ip(xi,x 2 )\ 2 . (55) 

The amplitude ([52]) leads to 

Pst = J dxdyf* n (x)W(x,y)f in (y), (56) 

where 

W{x,y) = Jdxidx 2 B*(xi,x 2 ;x) x 

[B(x 1 ,x 2 ;y) + B(x 2 ,x 1 -,y)]. 


/in(fc) = 


(49) 


(57) 
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The function (|53) ) gives 


W(x,y) = ±6(y-x)+^(3e r ^-e- r ' x -y\)e(-x)6(-y). 


The manifest effect of two photons to the system is 
realized by the initial wavepacket fi n (x) localized at x < 
0. The largest eigenvalue of W(x,y) gives the maximal 
probability of stimulated emission, and the corresponding 
eigenstate determines the shape of single photon wave- 
function fi n (x). 

Fortunately, the eigen-equation, i.e., the integral equa¬ 
tion 



dyW(x,y)f in (y ) 


A/in 0*0 


(58) 


can be solved exactly. In Eq. (1581) . the second order 
derivative to x leads to the differential equation 

dlfin(x) = ——^-T 2 fi n (x), (59) 

A 2 

where the solution is 

/ up2 \ \ V 4 / 7 

/in (4 = (wij ( 60 ) 


The eigenvalue A = 2/3 is obtained by the fact that fi n (x) 
is the solution of Eq. ([58]) . Finally, we conclude that the 
largest probability of stimulated emission in the two-level 
emitter is P s ™ ax = 2/3, and the corresponding incident 
wavepacket is f(x) = 2\/Te 2rx 0(—x). 

We notice that this result agrees with that in Ref. 0, 
where two types of the incident wavepackets, i.e., the 
Gaussian type and the Lorentzian shape are considered. 
By tuning the width of the wavepacket, the maximal 
probability 2/3 of the stimulated emission was found for 
the initial wavepacket with the Lorentzian shape. Here, 
we proved the exact result P™ ax = 2/3 by solving the 
integral Eq. analytically. 


B. JC system 

In this section, we consider the JC system as the sin¬ 
gle emitter coupled to waveguide photons, and apply the 
scattering theory to study the different scattering pro¬ 
cesses. The emitter Hamiltonian in the rotating frame 
reads 

H sys = A c a ] a + A e \e) (e| + g(a ] \g) (e| + H.c.), (61) 

where the detunings are A c = uj c — ko and A e = oj e — ko- 
For convenient, we focus on the resonant case A c = A e = 
0. The effective Hamiltonian H e ff = H sys — iTa)a follows 
from Eq. ([37]) . The external source term is 

Sj = "Vz r dt E w 

k,a=r,l 

+J k , a a\t)e- i ^ t ~ t d]. (62) 


For the single incident photon with the momentum k 
and the emitter initially in the ground state, the bound¬ 
ary condition is 7 i n = / and T i n = lim{j fc }_^ 0 d/dJ k , r in 
the asymptotic limit U —>• — oo. By the boundary con- 
dition 7out = I and F ont = lim|j >_j. 0 5/5J p ^p in the 
asymptotic limit tf —> oo, Eqs. (fMll and (pTHl lead to the 
reflection and transmission coefficients 

D X _ iTk X 

Rk6p ’~ k ~ (k + iT)k — g 2 ^ p ~ k ' 

Tk5p ’ k = (k + i r)/-p p>fc ' (63) 


For the two incident photons with momenta k\ and k 2 
to the emitter initially in the ground state, the boundary 
condition is 7 i n = / and T i n = lim^^o S 2 /SJ kl , r SJ k2 , r 
in the asymptotic limit ti —>• — oo. By the boundary 
condition 7 out = I and lim{j / j_ ) .o d 2 /SJ- Pl ^SJ- P2 ,i i n the 
asymptotic limit tf oo, Eqs. (flU) and (]3ll) lead to the 
8-matrix element 

c 7 7 _ p2 / + °° dt' 1 dt , 2 dt 1 dt2 ipi t\+ip 2 t , o-i kl t 1 -ik 2 t2 

* PlP2;fclfe2 “ (2 tt) 2 

x (Ta^a^a*(ti)al(t 2 )) , (64) 

of two reflected photons with momenta —pi and —p 2 - 
Here, the correlator for operators a(t) = e lH efft ae -iH eff t 
and a\t) = a ] e -iH eff t i ea q s to the 8-matrix ele¬ 

ment 


R—p\— p 2 ]k\k 2 R kl Rk 2 (S pikl 5 k2P2 + d P2kl 5 k2Pl ) 

r 2 g\E + iT)6 pl 

+P2,fel+fe2 

7r (E + iT)(E + 2iT) - 2g 2 
E(E + 2iT) - 4 g 2 

Yli=iA k i(ki + if) ~ g 2 ][Pi(Pi + iF) - g 2 }' 


(65) 


The Fourier transform of Eq. (l65l) leads to the wavefunc- 
tion 




■Mk 2 e 


ip\x\+ip 2 x 2 


P1P2 
iEx , 


_ e — c | R ki R k2 cos(kx) 

TV E s =±s(^-2A,)e i (f- A -) 3; 

" A+ - A_ (E + *r) (E + 2 iT) - 2 g 2 

: TFbIW+^wJ } (66) 


of two photons, where A s is the solution of k 2 +iVk — g 2 = 
0. As we discussed in Sec. Ill, under the weak driving 
light with frequency the quantum statistics of pho¬ 
tons emitting from the JC system can be described by the 
second order correlation function g^ 2 \x) = \rf(x c , x )\ 2 , 
where i/j(x c ,x) is the two photon wavefunction for the 
incident photons with frequency £q = &2 = uJd — • The 

above results agree with those in Ref. [22| • 

As a summary, in this section, we use the single emitter 
case as the example to show how our theory works in the 
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few photon scattering process. All the results agree with 
the previous studies, which justifies the validity of the 
developed scattering theory. 

V. VALIDITY OF MARKOV APPROXIMATION 

In this section, we use the simple example to study the 
condition of the Markov approximation in the array of 
two-level emitters. As shown in Secs. Ill and IV, the 
dynamics of the single emitter coupled to the waveguide 
with the linear dispersion can be exactly described by 
the effective Hamiltonian ([37]) and the generalized mas¬ 
ter equation dm. For the array of emitters, the dynam¬ 
ics is exactly characterized by the effective time-nonlocal 
action 5 e ff, where the effective Hamiltonian ([37]) and the 
generalized master equation ([26]) only describe the emit¬ 
ter evolution under the Markov approximation. 

In order to justify the validity of Markov approxima¬ 
tion, we compare the exact result given by the path in¬ 
tegral approach and the approximate result based on the 
quantum regression theorem. This comparison shows 
that the exact result and Markovian limit coincide when 
the bandwidth of the dynamics is sufficiently small com¬ 
pared to the distance between emitters. 

We focus on the single- and two- photon scattering 
processes, where the effective action 

S e « = f duYbt(u)[uI - Up {u)\ijbj(u) 

i-J 

~y j ( 6? ) 

j 

of emitters is given by the matrix with elements 
= -iTfSij ~ The driving 

action Sj = — f* f dtH ^ (t) is given by Eq. ([T9]) with 
Ok,da = The equations ([14]) 

and (131) lead to the exact results of the transmission 
spectrum and the second order correlation function by 
the effective action S e s and Sj. Under the Markov ap¬ 
proximation, Eqs. da and ([34]) give the approximate re¬ 
sults by the effective Hamiltonian Ho(w) ~ %o(fS) = Hq 1 
and O k ,a a ~ Oo,a a - 

A. Single photon processes 

For the single incident photon with the momentum k 
and the emitters initially in the ground state, the bound¬ 
ary condition is 7 i n = I and 7j n = lim S/5Jk, r in 
the asymptotic limit U — oc. It follows from Eqs. m 
and ([34]) that the reflection coefficient 

R k = -iT J2[G(k)] ij e i{k+ko ^ Xi+x ^ (68) 
ij 

of the photon with momentum p = — k is determined 
by the boundary condition 7 out = I and Tout = 



FIG. 3: (Color Online) (a) Single incident photon reflection 
spectra in a two-emitter system, where the distance between 
emitters is d ; (b) Excitation probability of second emitter 
|A(2,T)| 2 as a function of time T in a two-emitter system, 
when the second (first) emitter is initially prepared in the 
ground (excited) state. Here, T/ = 0 and T is taken as the 
unit. The Markovian results are also given by the solid (blue) 
curves. 

lirmr 8/8Jp,i m the asymptotic limit tf oo, where 
the Green function matrix G(k) = [k — Ho (fc)] -1 . In the 
Markovian limit, the reflection coefficient 

R f = —ir YyG u (k)]ije ik °( Xi+Xi \ (69) 
ij 

is determined by the approximate Green function 

G M {k) = lk-Hom~ 1 ' 

The difference of Eqs. ([68]) and ([69]) is that the phase 
factor e lkXi in Ho(k) is neglected in G u (k ), which re¬ 
sults in the condition kNd < 1 of the Markov approx¬ 
imation in the single photon scattering. In Fig. [3^, we 
show the reflection probabilities \Rk | and \R^ l \ of the 
scattering photon by two atoms with lattice spacing d, 
where kod = 2irn and n is an integer. The approximate 
result | R^jf | does not depend on d, which agrees with 
the exact result \Rk\ 2 very well for small kd <C 1. In 
the non-Markovian regime fed > 1, the dotted (black) 
curve (d = 1) shows that the reflection probability has 
the peaks localized around rioTr/d and no is an integer. 
The positions of these peaks are the resonant frequencies 
of eigenmodes in the “cavity” formed by the two emitters. 

The behavior of single excitation propagation in the 
emitter array can also be analyzed by Eqs. © and §1, 
where the exact result gives the condition of the Markov 
approximation. For the first emitter initially in the ex¬ 
cited state and the rest emitters in the ground state, the 
boundary conditions are Ti n = I and 7 i n = b\ at the ini¬ 
tial instant U = 0. The propagation of single excitation 
in the emitter array is described by the boundary condi¬ 
tion Tout = I and 7 ou t = bj at the final instant tj — T. 
The equations ([H]) and ([34]) lead to the amplitude 

^• r) 7£ (70) 

of detecting the excitation at the j-th emitter. The single 
photon emitted by the excitation of the first emitter is 
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described by the amplitude 


A a (x,T) = — iVT e~ 2 ' (7akoXj 0(cr a (x — Xj)) 

3 

x6(T - a a (x - Xj))A{j, T - a a (x - Xj)) (71) 


of detecting a single photon at the position x in the 
waveguide, where a denotes the right- and left- mov¬ 
ing modes, and the boundary condition is T 0u t = 
liuL{j fc }_^ 0 5/5Jk,(r,i) and Tout = I at the final instant 
tf = T. In the Markovian limit, the amplitudes become 


< 72 > 


and 


A™(x,T) = -iVT'}2e- ia ‘> koX i9(<T a (x - Xj)) 

3 

x6»(T - a a (x - Xj))A u (j,T- a a (x - xj)). (73) 


To investigate the non-Markov effects and justify the 
validity of the Markov approximation, we compare the 
exact results and the approximate results. It follows from 
the residue theorem that the poles £a and the correspond¬ 
ing residues Z\ of [1 /(c u — Ho(w))]ji determine the time 
evolution A(j, T) ~ Z\e~ l ^ xT . Similarly, A M (j, T) ~ 
T is given by the poles ^ and the corre¬ 
sponding residues Z^ of [l/(c j — 'Ho(^))\ji- Under the 
condition ^Nd <C 1, the phase factor e l ^^ Xi ~ x ^ ~ 1, 
and ^ and Z ^ are approximately the pole £\ and 
the corresponding residue Z\ of the exact propagator 
[1 /(c o — Ho(w))]ji- As a result, the condition of the 
Markov approximation is ^Nd <C 1. 

To understand this condition, we consider two emitters 
in the waveguide, where the first emitter is placed at the 
origin, and the effective action is determined by the 2x2 
matrix 


T-Lo(w) 


—iTf - iT -iTe i( < ko+ ^ d \ 
_ iTe i( k o+u;)d _ iTf _ iT J 


(74) 


It follows from Eqs. CUD and cm that the exact ampli¬ 
tudes are 

*4(1,T) = i e -( r /+r)q<7 + (T) + <7_(T)], 

*4(2, T) = l -e~^A^{C_(T) - C+{T)], (75) 

and 

A a (x, T ) = — iy/f [0(cr a x)6(T — a a x)A( 1, T — a a x ) 
-\-e~ teTak ° d d(a 0i (x — d))6(T — cr a (x — d )) 
xM(2,T - cr a x + cr a d)\ , (76) 

where 

OO 1 

C±(T) = -i[±r e iftod + (r /+ r)d (T - nd)] n 0{T - nd). 

n =0 

( 77 ) 
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FIG. 4: (Color Online) The solid (blue) and the dashed (red) 
curves show the wavefunctions of the right- and left- moving 
photon after the evolution time T, where the first (second) 
emitter is initially prepared in the excited (ground) state. 
Here, Tf = 0 and T is taken as the unit, and the green dots 
denote the two emitters. 


The amplitude A^(x,T) under the Markov approxima¬ 
tion has the same form of Eq. (ESI, where A(j,T) is 
approximated by 

A M (1,T) = cosh(rTe iA;o<i )e“ (r ^ +r)T , 

A u (2, T) = - sinh(rTe ifeod )e“ (r ^ +r)r . (78) 

The agreement between Eq. (ESI and the exact result 
([75]) is verified in the Markovian limit (Tf + T)d <C 1, 
where C±(T) ~ 

The time evolution of the two emitters exhibits two 
non-Markovian effects. The first is the retardation ef¬ 
fect. The amplitude M m ( 2, T) shows that once T > 0 the 
second emitter has the probability in the excited state. 
However, this is an artificial effect of the instantaneous 
effective Hamiltonian H e q under the Markov approxima¬ 
tion. In fact, the single photon wavepacket emitted by 
the first emitter takes the time T = d to arrive at the sec¬ 
ond emitter. Hence, the second emitter has to stay at the 
ground state for T < d, i.e., A( 2, T < d) = 0. This retar¬ 
dation effect is fully characterized by the exact result dm 
and C±(T). We show the probability \A(2,T)\ 2 in Fig. 
Eb for Tf = 0 and kod = 27m, where the dotted (black) 
curve for d = 1 displays the retardation effect explicitly. 
In the small d limit, i.e., Td/c <C 1, the Markovian result 
agrees with the exact result very well, as shown by the 
solid (blue) and dashed (red) curves in Fig. [3 Jd, where 
the speed of light c is taken to be the unit. 

The second effect is the formation of the entangle state 
(| eg) — \ge))/y/2 in the limit T —>• oo. For the vanishing 
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FIG. 5: (Color Online) (a)-(b) The exact and Markovian T- 
matrix elements for E = 1; (c)-(d) The second order corre¬ 
lation functions of out-going photons: In (c) r = 0, in (d) 
E — 0, and Markovian results are shown by the solid (blue) 
curves. Here, Yf = 0 and T is taken as the unit. 


Yf = 0, both the exact and Markovian results show that 
in the steady state T —>• oc the emitters have the proba¬ 
bility P e to form the entangle state. Under the Markov 
approximation, the probability P e = 1/2 is determined 
by M m (1,oo) = —M m (2,oo) = 1/2. The exact result 
M(l,oc) = — *4.(2, oo) = 1/(2+ 2 Yd) gives the probability 
P e = 1/[2(1 + Yd) 2 ] to form the entangle state. In the 
limit Yd ^ 1 , the Markov approximation works perfectly, 
i.e., A(j, oo) ~ *4 M (j,oo) = 1/2. 

During the formation of the entangled state, the dy¬ 
namics of the waveguide photon is described by the am¬ 
plitudes A r ,i(x,T). In Fig. m we show the propagation 
of the right- and left- moving wavepackets by A r ,i(x,T) 
for the distance d — 1. In the steady state T ^ oc, the 
entangle state is established, where the standing wave in 
the regime 0 < x < d forms to mediate the interaction of 
two emitters. 


B. Two photon processes 


In this section, we study the scattering process of two 
incident photons with momenta k\ and &2 • By comparing 
the exact result and the approximate result, we investi¬ 
gate the condition of the Markov approximation in the 
two-photon processes. 

The initial boundary condition is y in = I and P- in = 
lim{j fc T.^ 0 b‘ 2 /fiJk 1 ,rfiJk 2 ,r in the asymptotic limit U —>• 
—oo. It follows from Eqs. ([Hj) and (15311 that the S- 


matrix 


^PlP2,klk2 i4li4 2 (Wp2*2 "t" ^P2kl 3pik2 ) (79) 

F 2 

— i ^Pi+P2,&i+&2 ^ ^ Gi{pii P2)Tij{P)Gj {kl, kyd) 

7r z ' 

*j 

of two reflected photons with momenta —p\ and — P 2 are 
determined by the boundary condition 7 out = I and 
Fo Ut = limijj. 5 2 /5J- pit i5J- P2 j in the asymptotic 
limit tf —> oo. Here, 

G i (p 1 ,p 2 ) = 'E G hi(P i)G i2i (p 2 )e i ^ +ko)x ^e i ^ +k ^, 
*1*2 

(80) 

and the T-matrix T(E) = —II 1 (E) is given by the 
Dyson Expansion with the bubble 

I lij(E) = i J ^-Gij(w)Gij(E - to) (81) 

and E = k\ + 

Under the Markov approximation, S PlP2i k 1 k 2 is given 
by the same form of Eq. o, where Rk i ~ R™., 

Gi(pi,P2) ~EC(Pi) G “fe)e ito(lil+li2) (82) 


and the bubble (1811) are determined by the approximate 
Green function G m (ci/) ~ [uj — Ho(O)] -1 . The exact result 
shows that the Markov approximation is valid in the limit 
kiNd , piNd , and YNd <C 1. 

In order to understand the condition, we investigate 
the two-photon scattering processes explicitly by consid¬ 
ering the photon scattering by two emitters, where kod = 
27m. We first compare the exact T-matrix T(E) with 
the T-matrix T M (E) under the Markov approximation, 
where the elements 11^ (T) = II^ 2 (T) = II^(T) + F[jf(T) 
and U™(E) = n ™{E) = II™ (E) - II™ (E) are given by 


n “(^) = J E 

cr=±l 


1 

E + 2iYf + 2iT(l + ae ikod ) ’ 


U™(E) 


1 1 

2 E T 2iiYf T 2iT 


(83) 


The absolute values of the exact T-matrix elements are 
shown in Figs. [5^ and [5 Jd for kod = 27m, which illus¬ 
trate the perfect agreement between the Markovian re¬ 
sult T™(E) and the exact result in the limit Ed , Yd <C 1. 

The Fourier transformation of S'-matrix leads to the 
wavefunction 


'ip(x c ,x) = I J dp±dp 2 S PlP2 , kl k 2 e lpiXl+lp2X2 (84) 

of two reflected photons, where the center of mass co¬ 
ordinate x c = (xi + # 2)/2 and the relative coordinate 
x = x\ — X 2 - For two incident photons with the same 
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momentum £q = k 2 = fc, the photon statistics is charac¬ 
terized by the second order correlation function 


9 {2) (t) 


I^Qect)! 2 

I^l 4 


(85) 


In Figs. [5b and [5ji, we show g( 2 )(0) as the function 
of E = 2k and g^ 2 \r) for the resonant frequency k = 
0, where the Markovian results perfectly agree with the 
exact result in the small kd and Yd limit. When d is 
increasing, e.g., Yd = 0.1 and 1, the Markovian result 
deviates from the exact one. 

As the summary in this section, we use the exact result 
to study the non-Markovian effects and examine the con¬ 
dition of the Markov approximation. We conclude that 
in the limit YNd , kNd <C 1, the Markov approximation 
works perfectly. For the large YNd and kNd, some non- 
Markov effects emerge, e.g., the retardation effect. In 
the following sections, we assume the system length Nd 
is small enough and the Markov approximation is always 
valid. 


VI. ENTANGLED PHOTON PAIRS BY SINGLE 
EMITTER 


In this section, we use the scattering theory to study 
the generation of entangled photons in the scattering pro¬ 
cess. In order to realize the deterministic generation of 
reflected photons, the two-level emitter is placed at the 
left hand side of a perfect mirror, i.e., in the half-end 
waveguide, as shown in Fig. [lb. Here, the mirror is put 
at the origin, and the position of the emitter is xo < 0. 

It follows from Eq. ([35]) that the effective action is 


S eff - S sys + i j dtTe\t)e{t)+i j dtY b b\t)b{t) 
-\-iy/YY b J dtb\t)e(t + xo)e~ lk ° x ° 

+i^/YY b f dte\t)b(t + xo)e~ lkoX °, (86) 


where the emitter action S sys is determined by the Hamil¬ 
tonian H sys = ^emitter• By the creation (annihilation) 
operator e t (e), the hardcore boson is introduced to de¬ 
scribe the two-level emitter, where the energy level spac¬ 
ing is uj e = fco and ^emitter = Uoe^e^ee/2 in the limit 
Uo oo. The external source term Sj = — f* f dtH&(t) 
is given by Eqs. m and (J6]) with the jump operator 
O = e. 


condition 7 out = / and E out = lim { Jk y^ 0 6/6J Pi i in the 
asymptotic limit tf oo, Eq. (fl^l) and (|34|) result in the 
reflection coefficient 

RkS P ,- k = -iS Pt - k [TG ee (k)e 2 ^ k+k ^ X0 +T b G bb (k) 

+20T ~ b G eb (k)e i{k+ko)x °}, (87) 

where the Green functions G ss ’(k ) = —i f dte lkt (s{t)s' Xy ') 
with s, s' = e, b. The effective action determines the sin¬ 
gle particle Green functions 


G ee (k) G eb (k ) \ 

G b e(k) G bb (k) ) 
k + iY 


i^YYle-^+b^ k 
which give the reflection coefficient 


iY b 


Rk3p, — k — 


k — iY iYe~ 2i ( k + ko ^ Xo \ 


r 

k + iY - iYe 2 d k o+ k )\xo\ °v~ k 


,( 88 ) 


(89) 


and \Rk\ = 1 in the limit Y b oo. The phase shift 
arg (Rk) is shown in Figs. [6b for d = 10 -4 , where 6q = 

fco \xo\- 

For two incident photons with momenta k\ and k 2 and 
the emitter initially in the ground state, the boundary 
condition is 7 i n = I and T \ n = lim{ Jk y^o S 2 /SJk 1 , r fiJk 2 ,r 
in the asymptotic limit U — oo. By the boundary con¬ 
dition 7out = I and J ~;>ut = lim^j fc |^.o^ I^J—pi,i^J—p 2 ,i 
in the asymptotic limit U — oo, the S'-matrix element 


^—Pi—P2\k\k 2 — RkMSp^Sp^ H - 3pik 2 fipik 2 ) 


.16 


i Y 5 Pl+P2: k 1 +k 2 T{E) x 


(90) 


n 


sin [(kj + fc 0 ) \xp\\ sin[(pj + fc 0 ) \x Q \\ _ 

'—i^2 + ^F — iYe‘ 2 d k °+ ki ')\ Xo \]\pi + iY — ire^^o+pdkol] 

of two reflected photons with momenta p\ and P 2 is given 
by Eqs. m and (O in the asymptotic limit tf oo. 
Here, the T-matrix element 


T(E) = - w 

v ’ ur 1 


l 


o -"(/•-") 


is determined by the vacuum bubble 
duo 


n (E) = i j 


2i r 


G ee {ui)G ee (E — uj) 


(91) 


(92) 


and E = ki + k 2 = pi + P 2 - 


B. Entangled photon pairs 


A. Single- and two- photon scattering 

For the single incident photon with the momentum k 
and the emitter initially in the ground state, the bound¬ 
ary condition is 7 i n = I and E[ n = lim ^y^S/SJ^r 
in the asymptotic limit U — oo. By the boundary 


For the state E fclfea f iki) f {k 2 )r\ i r\^ |0) /\/2 of two in- 
dependent incident photons, the ^-matrix dm leads to 
the asymptotic state 

|VWt) = “y= V ; out(Pl:P2)^- Pl ^-p 2 |0) (93) 

V Z P1P2 






15 



FIG. 6: (Color Online) (a) The phase shift of the single pho¬ 
ton; (b) The Von-Neumann entropy of the reflective photons. 
Here, d = 10 —4 , and T is taken as the unit. 


of two reflected photons, where the wavefunction is 

^out(pi,P2) = f(pi)f( P 2)R Pl R P2 - 16iT 2 T(E)F 2 (E) 
tt sin [{pi + fc 0 ) |xq]| 

-I-L p • 4- iT — ire 2< ( fc o+Pi)l*o| 5 

i= 1,2 11 


and the integral 


F ( E )= f — TT /(f + crfc ) sin[(f +crfc + fco) |a?o|] 

2 J ^ I +crk+iT-iTe 2i ^ +crk+ko '>^' 


(95) 

The Von-Neumann entropy SVn = — tr(p\\npi) is 
given by the single-photon reduced density matrix pi = 
^ 2 (|Vwt) (t^out |) i where the degree of freedom for the 
other photon is traced out. The Von-Neumann entropy 
Svn = — l n can be obtained by the singular 

value decomposition of V’outbi,^) = J 2 \ 9 \(Pi)^ 9 \(P 2 )- 
Here, the singluar values A measure the entanglement of 
two photons in different modes g\ and g\. For the single 
photon wavepacket 



1 

k + iy’ 


(96) 


In the rotating frame, the Hamiltonian © becomes 
H sys = ^[A e e|ei + + (fiefs, + H.c.)] 

i 

T-f^HC + 2 y! Uijs\s^SjSi. (97) 

ij 

where A e = uj e — &o, A s = uj s — ko + ujd, and we focus 
on the two photon resonance case A s = 0. The waveg¬ 
uide photons couple to the N atoms collectively through 
the operator O/ Cj ± = y/T e~^ k±k °^ Xi e*/ y/L , and the 

free space modes couple to the atom operator locally, 
which induces the decay of the excited state e. The ef¬ 
fective action 


<S e ff = J dt(y' j [el(t)(idt+iTf)e i (t) 

+4 (t)id t Si(t)] - H sys } + Sre, (98) 


of EIT atoms is given by 

S re = iT f dt e*(t)ej(t — \xj — Xj\)e lk °^ Xi ~ Xj ^. (99) 
ij 

The source term Sj = — f* f dtH&(t) is determined by 
Eq. ([T9 j) . Here, the lattice spacing d satisfies kod = 
(2 n + 1/2)7t. 

Based on the general result da and fl, we first 
study the single- and two- photon scattering, and show 
the transmission spectrum and the second order correla¬ 
tion function of out-going photons. In the second part, 
we investigate the transient processes and show how the 
wavepackets of single and two incident photons transfer 
to the atom excitations, propagate in the Rydberg-EIT 
atom array, and finally emit back to the waveguide. For 
the single incident photon, we show the free propaga¬ 
tion of the dark polariton. For two incident photons, we 
show the counter-propagation and co-propagation of two 
polaritons. 


we show Sy n in Fig. [ 6)3 for d = 10 4 , which displays the 
generation of entangled photons by the two-level emitter. 


VII. RYDBERG-EIT SYSTEM 

In this section, we consider the application of the scat¬ 
tering theory in the photon transmission in the EIT 
atoms coupled to the Rydberg level. We shall show that 
the developed scattering theory is an efficient approach to 
the photon transmission in the array of interacting emit¬ 
ters with complicated structures. Here, we highlight the 
interplay between the EIT phenomenon and the Rydberg 
interaction, and show a rich variety of quantum statistics 
of the scattering photons and polariton exitations. 


A. Single and two photon scatterings 

For the single incident photon with the momentum k 
and the atoms initially in the ground state, the bound¬ 
ary condition is 7 i n = / and = lim{ fi/$Jk,r in 
the asymptotic limit ti —»■ — 00 . By the boundary con¬ 
dition 7 out = I and JVut = lim {j k }^oS/SJ p ,( r ,i) in the 
asymptotic limit tf 00 , Eqs. (|m|) and (IM1) lead to the 
reflection and transmission coefficients 

RkS P ,- k = ~iT6 P! - k Y^[Go(k)]tJe^ k+k ^ Xi+x R,( 100) 

ij 

T k S p , k = 5 pk {l-iYY}G a mt^~ i(k+ko){xi ~ Xi) }, 

ij 

of the photon with momentum p. 
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FIG. 7: (Color Online) The single photon transmission spec¬ 


tra, where T = l, A e = 0 , H = l, and the atom number is 
20. (a) Tf = 0 and d = 10“ 4 ; (b) T f = 1 and d = 10 -4 ; (c) 
T f = 0 and d = 10 —2 ; (d) T f = 1 and d = 10 -2 . 


The free Green function 
GoM = 


1 


uj - 

of atoms is determined by the block form 


( 101 ) 


nn(u) & f(Ae - iTf )6 ij - ire^o+-)|x i -x,| n5ij \ 

0[ ’ \ nSij a Aj J 

( 102 ) 

in the basis |e^) and |s$), where each matrix element is 
the TV-dimensional matrix in the coordinate basis. In 
the notation [Go(k)]fj T , cr(a') = e,s denotes the different 
“spin” blocks and i, j denote the coordinate in the block. 
In the Markov limit uod 1, T-Lq{(jS) is approximated by 
the frequency independent effective Hamiltonian 


njM 

n 0 


(A e - iT f )Sij - iTe ik ^ Xi ~ x ^ nS tj \ ( . 

a Aj J ’ 1 J 


In Fig. El the exact transmission probability is com¬ 
pared with that under the Markov approximation. As 
shown in Figs. EH and Eb, for the small lattice spac¬ 
ing d = 10 -4 , the Markov approximation works very 
well, which leads to the same result as that from Eq. 
(1102D . When the lattice spacing d = 10 -2 is larger, Fig. 
Et shows the difference between the results under the 
Markov approximation TO and the exact result (IT 02 lh 
In the realistic case, the small length of the array jus¬ 
tifies the validity of the Markov approximation. Hence¬ 
forth, we focus on the Markov limit. In Figs. Et> and 
EH, we have taken into account the decay of the excited 
state | ei) to the free space, i.e., Tf ^ 0. The single pho¬ 
ton transmission shows the EIT nature of atoms, where 
the total transmission appears at the resonant frequency 


k = 0 in the EIT window. The reason for the total trans¬ 
mission is that the resonant photon transforms to the 
free-propagating dark polariton excitation, which is only 
the superposition of states r\, \gi) and | Si). As a result, 
the free space decay T f of the state |e$) does not affect 
the total transmission of the resonant dark polariton with 

fc = 0. 

For the two incident photons with momenta k\ and 
&2 and the atoms initially in the ground state, the 
initial boundary condition is 7 i n = / and T \ n = 

lim{j fc }_^o d 2 /dJk 1 ,rdJk 2 ,r in the asymptotic limit U 
—oo. By the boundary condition of the final state 
7out = I and J'out = lim{j fc }_> 0 5 2 /6J Pur 5J P2:r in the 
asymptotic limit and tf oo, Eqs. (ITU) and (IM1) lead 
to the two photon ^-matrix 


^PlP2]kik2 — Tk 1 Tk 2 (dp 1 k 1 fip2k2 fipik2^P2kl) 




(27r) 2 . ^ . 

U« 2 ,jlj 2 


O iko (xj t +Xj 2 —Xi t —Xi 2 ) 


(104) 


UP2-MM) 


for two transmitted photons with momenta pi and p 2 - 
Here, the Green function 


G i e £- jlj2 (P uPXhM 


= J dt[dt , 2 dt 1 dt 2 e ipi ^ +ip2t2 - ikltl - ik2t2 

x (t^e) (t 2 )) (105) 


is the Fourier transformation of the four-point connected 
Green function. 

The Dyson expansion in terms of the two-body inter- 
action Hue + J2ij Uijs\sjSjSi/2 leads to 


(o) r 2 

SpiP2]k 1 k 2 = ^Pip 2 'Mk 2 ^9 ^Pi+P2M+k 2 x (106) 


E 

iji'j' 

<J\(7 i<J2<^2 

+ (Pl ^P2), 


27r 

[^(pi,P2)}T [ inm-ti' 2 *' 2 Mh,k2)]?r 


where we define the independent scattering part 
S Pi P 2 -Mk 2 = T kiTk 2 d Pl k l 5p 2 k 2 , and the function 


Mki,k 2 W; = J2e iko{x “ +x ^[G 

j 132 

( 10 7) 

The T-matrix, depicted by the ladder diagram in Fig. 
Eb, satisfies the Lippmann-Schwinger equation 


u <J 1 <T2 ^Crjcr^ 


E [Tmz&r 


(108) 


+u aia 




where the vacuum bubble is 

mm&f = i j^[GoMwiGoiE-ur;/. ( 109 ) 
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In the matrix form, the Lippmann-Schwinger equation is 
formally solved as 

T < £ > = u-.-rw < 110 > 

where the vacuum bubble IIo(^) = (E — is given 

by H .2 = and the interaction matrix 

U has the diagonal element Uf? = U?? = U£? = Uo ^ oc 
and Uff = Uij in the basis {| eiej) , | eiSj) , | s&j) , 

The wavefunction 

Ip{x c ,x) = J ^Els piP 2 ., klk 2 e i ^+ i P— 

AEx c 

= -^ r {2T f+fc T f _ fc cos(M-^ ^ [F(x)]?f 

iji'j' 

<j<j' nn' 

x [T(E)]™.;f [w (| + k, | - k)]?:,}, (Ill) 

of two transmitted photons is the Fourier transform of 
Sp 1 p 2 -,kik 2 ■> where 


icrcr 

\ij 


[F(x)} 

x*( icr )x;*(i cr ')[ e * (f 

+ (x —>> —x). 


_ p *^0^c,12 _ 

i E p-Xi(ne)Xi'fee) 

“ A - £1 - £v 


1*2 

hf -ez')® 


0(x) + e 


-i(fsi)xn(_ 


9(-x)} 


( 112 ) 


is determined by x c? i 2 = +%i 2 and the eigenstates |%/) 

(| xi)) of Hq 1 with the corresponding eigenenergies 

si (e z *). Here, (ia\xi) = Xi(i°) and I Xi) = Xi( i(T )- 
By the wavefunction (jTTTj) , we show the nor¬ 
malized second order correlation function g^ 2 \x) = 

2 

7T0(a; c , x)/Tfli 2 of out-going photons for the two in¬ 
cident photons with the same momentum k\ = &2 = E /2 
in Figs. Eb and [8fc* For the resonant case E = 0, Fig. 
Eb shows the second order correlation functions for the 
uniform interaction Uff = C, the van-der Waals inter¬ 


action Uff = C§/\i — j| 6 , and the dipolar interaction 

Uff = Cs/ \i — j| 3 , which exhibit the anti-bunching be¬ 
havior of out-going photons. 

It can be understood in the following way. As shown in 
the left panel of Fig. Etb if two Rydberg excitations are 
close to each other, the Rydberg state is shifted by the 
strong Rydberg interaction, such that the classical light is 
off-resonance with respect to the transition between |e*) 
and | Si). As a result, the photon with k = 0 is resonant 
with the transition \g/) |e$) and reflected. For two 

transmitted photons, to maintain that the frequency k = 
0 of the photon is in the EIT transmission window, they 
repulse each other and show the anti-bunching behavior 
such that the Rydberg state is not shifted. 

For the incident photons with frequency E/2 larger 
than some critical value, the transmitted photons can 
also display bunching behavior, as shown in Fig. Et- The 
mechanism of the generation of bunched photons is illus¬ 
trated in the right panel of Fig. Eb- For two Rydberg 





FIG. 8: (Color Online) Bunching and anti-bunching behaviors 
of transmitted photons, where r = 1, Tf = 1, A e = 0, H = 1, 
kod = 7r/2, Uo = 10 8 , and the atom number is 20: (a) The 
Feynman diagram for the T-matrix; (b) The second order 
correlation functions g^ 2 \x) for C — C 3 = Cq = 1, E = 0, and 
k — 0; (c) The second order correlation functions In 3 ^( 0 ) 
for the uniform case C — 0.46, where the relative momentum 
k — 0, and the inset shows g^ 2 \x) for the frequency E/2 — 
0.2 of each photon; (d) The schematic for the generation of 
the bunched and anti-bunched photons, where the incident 
photons have different frequencies, s' is the shifted energy 
level due to the Rydberg interaction. 


excitations close to each other, the Rydberg energy levels 
|are shifted. The photons with finite momentum E/2 
realize the two photon resonance with the shifted energy 
level |Si), and can be transmitted. In order to achieve 
the transmission of photons by the shifted Rydberg lev¬ 
els |Si), the two Rydberg excitations prefer to stay next 
to each other, which induces the bunching behavior of 
the transmitted photons. 

As discussed in Sec. Ill, the second order correlation 
function g^ 2 \x) obtained by the scattering theory can 
characterize the quantum statistics of photons emitted 
by the Rydberg-EIT atoms under the weak driving field. 
In Ref. (36j| , the result from the scattering theory and the 
numerical solution of the master equation for the effective 
spin model under the weak driving light are compared, 
where the two results agree with each other perfectly. 


B. Propagations of single and two excitations 

In this section, we investigate how the single and two 
incident photons transform to the excitations of the Rdy- 
berg atoms, and the propagation of excitations. 

For the single incident photon with the wavepacket 
f(k) and the atoms initially in the ground state, 
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(b) 

A(T) = free 

^ propagation 


+ 



polariton 

scattering 


FIG. 9: (Color Online) (a) Probability of Rydberg states for 
the single incident photon wavepacket, where r = 1 , Tf = 0 , 
A = 0, Q = 0.1, d — 10 -4 , 7 — 0.01, and the atom num¬ 
ber is 20. (b) The Feynman diagram for the two excitation 
propagation. 



the boundary condition is 7 i n = / and T \ n = 

S t /( r ,i)(fc)lim{4}^oV w t,('-,i) at the instant U = 0, 

where the wavepackets of the right- and left- moving pho¬ 
tons are 


/W = /r(*) = J^rr, (113) 

V 7T k + t'y 

and 

/ /y p ikx N 

- 1 —, (114) 

7T k — 27 

with the width 1 / 7 . By the final boundary condition 
T ,out = 1 and 7 ou t = i±i = e^Si at the instant tf = T, 
Eqs. (fl4|) and ([34]) lead to the amplitude 


A^iT) = iVoVArf^e^^eiTa - <r a Xj) X (115) 
lj 

U e ) r -7(T a -q a a:,-) e~ i£l ( Tcx ~ (TaX j)] 

Si + 27 

of the excitation g = e, s at the position i, where T a = 
T — z a and 27 = 0, zi = xy. For this situation, the dark 
polariton forms, where the probability of the occupation 
in the excited state \ei) is quite small ~ 10 -4 . In Fig. 
[9l we show the dark polariton propagation for the single 
right moving incident photon. Tn order to show the slow 
propagation of the dark polariton, we choose the small 
ratio Q/T = 0.1. 

For the wavepacket of two incident photons and the 
atoms initially in the ground state, we consider both the 
co-propagation and counter- propagation cases. For the 
co- propagation, the initial boundary condition at the 
instant U = 0 is 7 i n = I and 


^in 


£ fr(kl)f r (k 2 ) 

k\k,2 


lim 


£ 2 


{-h;}—>-0 SJ kur SJ k2 


(116) 


while for the counter- propagation, the initial boundary 
condition at the instant ti = 0 is 7 i n = I and 


^in 


£/r(*l)/l(k) 

kik2 


{Jfc}—>*0 SJk 1:r SJk 2 ,l 


(117) 


FIG. 10: (Color Online) Probabilities in the state |s) for two 
counter- and co-propagating polaritons with the dipolar case 
C 3 — 1, where F = 1 , Tf = 0, A e = 0, d = 10 —4 , and 
the atom number is 20 . (a)-(c) show the probabilities at the 
instants T = 600, 1200, and 1800 for the counter-propagation 
case: Tl = 0.1 and 7 = 0 . 01 ; (d)-(f) show the probabilities at 
the instants T = 12, 20, and 28 for the co-propagation case: 
Q = 1 and 7 = 0 . 1 . 


By the boundary condition T 0 ut = 1 and 7 ou t = 
Mi,ii/i 2 ,i 2 (Mi ,2 = e, s) at the instant tf = T, Eqs. (fT41) 
and ([34]) result in the amplitude 


= V m „{a£1(T)a£L(T) 

+2~iTa a £ e ik ° x n +iff “ feo ^ e(T - x h ) x 
3lh 

9(T -z a - a a Xj 2 ) (118) 

E [_1_ y 

}w 1 +W 2 -H2 iili2 ’ ii ij 

13 ,00 


[G(u>i)]ij [G(u> 2 )]? 


a e 
332 ‘ 


0 -iuj 1 (T-x jl ) p -iu 2 (T-z a -(T a Xj 2 ) 


0 J 1 + 27 


Cl>2 + ^7 




of two excitations Mi and M 2 at the positions i\ and 12 for 
the co-propagation a = r and the counter-propagation 
a = l case, where the operator Vi llx ^i 2 ^ 2 symmetrizes 
the wavefunction under the interchange Hmi <—> £2 M 2 - 

In the first row of Fig. we show the propagation of 
two dark polaritons for two counter-propagating incident 
photons. Here, we choose fi/r = 0.1 to show the collision 
of two slow polaritons. When the two excitations ap¬ 
proach with each other, the states |s$) are shifted by the 
Rydberg interaction, which results in the off-resonance 
with the transition from \ei) to | Si) induced by the clas¬ 
sical light. As a result, the photon is resonant with the 
excited state |e$) and reflected. After the collision of two 
excitations, they propagate away from each other, and 
the interaction gradually vanishes, which results in the 
re-formation of two free propagating dark polaritons. 

In the second row of Fig. \W\ the propagation of two 
dark polaritons for the co-propagation incident photons 
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is shown. When the first photon transforms into the dark 
polariton in the atom array, it blocks the transmission of 
the second photon. This blockade occurs over a char¬ 
acteristic distance r& ^ (CsT/D 2 ) 1 / 3 . For an incoming 
two-photon wavepacket whose size is larger than r&, the 
blockade manifests itself as a suppression of the proba¬ 
bility of two photons to overlap with each other as they 
propagate through the medium (see Fig. 10 d-f). During 
the propagation in the Rydberg-EIT atoms, the two po- 
laritons keep away from each other such that the energy 
level | are not shifted, which results in the transmission 
of two anti-bunched dark polaritons. 


the dissipative many-body systems. 

Acknowledgments 

This project has been supported by the EU project 
SIQS. DEC acknowledges support from Fundacio Privada 
Cellex Barcelona, ERC Starting Grant FOQAL, and the 
Ramon y Cajal program. Tao Shi thanks the useful dis¬ 
cussions with Y. Chang, A. G. Tudela, C. N. Benlloch, V. 
Paulisch, C. S. Munoz, T. Caneva, and M. T. Manzoni. 


Appendix A: Generalized master equation by path 
VIII. CONCLUSION integral 


We summarize our results in this section. We devel¬ 
oped the scattering theory to investigate propagation of 
photons through an array of quantum emitters using the 
path integral approach. The exact transition amplitude 
for arbitrary initial and final states is obtained to de¬ 
scribe the quantum statistics of scattering photons and 
the dynamics of emitters in the transient process. The ex¬ 
act result justifies the correctness of the Markov approx¬ 
imation for the single emitter coupled to the waveguide 
photons with linear dispersion. The exact and Marko¬ 
vian results coincide when the bandwidth of the dynam¬ 
ics is sufficiently small compared to the distance between 
emitters. Here, the generalized master equation for the 
few-photon scattering process is obtained to describe the 
transient dynamics of emitters. The generalized master 
equation establishes the relation between two equivalent 
systems, i.e., few photon scattering by the emitters and 
the emitters under weak driving light. 

For the single emitter case, two paradigmatic exam¬ 
ples, i.e., the two-level emitter and the JC system, are 
used to show the correctness of our theory by compari¬ 
son with the well-known results. For an array of emitters, 
the validity of the Markov approximation is examined in 
the system with two-level emitters coupled to waveguide 
photons. Here, the dynamical evolution of emitters show 
some non-Markovian effects, i.e., the retardation effect. 

The generation of entangled photons by the single 
emitter in front of the mirror is also analyzed by the 
exact result from our theory. Finally, the photon trans¬ 
mission in an array of EIT atoms coupled to the Rydberg 
level is investigated by the scattering theory. We high¬ 
light the interplay between the EIT phenomenon and the 
Rydberg interaction, and show how this results in the 
bunching and anti-bunching behaviors of the scattering 
photons and the dark polaritons propagating in the array. 

The scattering theory also provides the way to explore 
some non-Markovian effects of emitters coupled to the 
waveguide with nonlinear spectrum, where the multi¬ 
photon bound states (HU, EBj may form. The general 
result of the transition amplitude and the generalized 
master equation enable us to study the photon transmis¬ 
sion in more complicated quantum optics systems and 


In this Appendix, we derive the generalized master 
equation ([26]) by the path integral approach. To sim¬ 
plify the notation, we neglect the term Hr f describing 
the independent decay of each emitter to the free space, 
and the derivation including the free space mode follows 
the same procedure. In the end, we show the result by 
considering the effect of the free space decay. 

Here, we need to introduce the emitter field explicitly. 
In quantum optics systems, the emitter operators are 
usually the annihilation (creation) operators of bosonic 
modes and the ladder operators. If one use the hard¬ 
core boson to describe the emitters, all the operators 
are bosonic operators (6/), where l denotes the different 
modes of the emitters. In the coherent basis, the re¬ 
duced density operator p s (T) = Tr^\\e~ lHT ppd)e lHT } 
becomes p s (T) = ^n J -^ r in,+ e ^ fe,Q! Jk ’ a ~ Jk ’ a ' + pj(T): 

Pj(T) = Jdp(Pi r + iOU tiPi i + iOU t)dP'(Pi,—,outiPi i _ jOU .t) 

e~ A,-,out; T) 

X |{A,+,out}> ({A,-,out} | , (Al) 


where a = r, l denotes the right- and left- moving photon 
modes, 


? in,± = lim y: ^in({rcfcq})TT 

\ Jhoc J 


S nka 


f T^ka .} 


ka 




k,a,± 


(A2) 

the measure for the unnormalized coherent state \{Pi}) 
is dniPuPf) = n*( e d 2 Pi/tt), and the element 

Pj(Pt+,out,Pl,-,ouf,T) = j dMA, ± , in) tf, ± , in ) X 

^(A r +,out5 A, — ,out? A, + , in? A 5 -,irb ^)p0 (A, + ,in 5 A,-,in) 
is given by the propagator 


^AA, + ,out? A, — ,out5 A,+,in 5 A, — ,in5 T) 

= J dp(Jk,a.i Jk,a) (A 3 ) 

X ({J K a}\ <{A, + ,out}l e~ iHT |{A, + ,in}> | {Jk,a,+ }) 

X <{J fc|a ,-}| ({A, — ,in}| e iHT |{A,-,out}> |{Jfc,«}) 
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in the closed time path. Here, we define 

A)(A,+,in’ Pl,~ ,in) = ({A, + ,in}| Psys(0) |{A, —,in})- 

The propagator (|A3| ) is obtained by the saddle point 
method shown in Sec. IIIC, which is 

Z j(Pl,+,onU A,-,out; A,+,in, A*-,in! O 


e* 56 x 


J d/l(Jk,a, Jk,ct) 


ex 


p{ E {Jt a Jk,a, + e- ia - kT + JU.J k , a e^ kT ]} 


k,a=r,l 


J D [system] ( 


.iSeff ,+ —i-S'eff _ + iSj,+ — ZSJ _ 


(A4) 


Here, iS b = Ei [A*+,outA(0 + A*CO A,-, out] is the 
boundary term for the emitters, and S e ^^± are the effec¬ 
tive actions in the forward and backward time-evolution 
paths, which are given by substituting the fields /3^± for 
the emitter fields in the effective action S e fi . The external 
source term is 


Sj,± = - f dtY}JZ, a Ok,„ a , ±(t)e-~“ fc(r - t) 

J ° k,a 


+J k . a ,±0* k ^ ± (t)e— kt } 


(A5) 


Finally, the Gaussian integral over J \ CjQ ,, J £ a leads to the 
propagator 

^(A,+, out’ A,-,out; /?/, + ,in J A, — , ini 

= J J D[system]e^' s,d ’ +_ ' s,d ’- + ' s,jump ^ ) 

exp{E[A: +) outA(T) + A*(OA,-, out]}, (A6) 

l 

where 


AtJoikt 


S d>+ = S eff>+ - / dty[4,_O t ,a 0 , + (f)e l 2 3 4 5 6 7 

Jo fc, Q 

+J fc)a , + 0^ )+ (t)e-^ fct ], (A7) 


Acr a kt 


5 d> _ = s e y_- /U*E[4Y-°^,-(*K' 

- 70 fc.a 

+^,a,+Ofc ) .„,-We- iCTafct ], (AS) 

and the jump term 

s iump = -i j dtdt'^e^-^Ok^A^oi^An 

k,a 

(A9) 

Under the Markov approximation, the actions read 

ftf 

xO* ± {t)O h± {t)e ik ° (A10) 


S e S,± = S'sys.i +i [ dt^2 

Jti 


and 


Sjump 


rT 

/ ^ \ZrjTj cos[fcp(^ — ; 

Jo 


xO,,+ (t)0*_(t). 


(All) 


By taking the time derivative of Eq. CEQ> and using Eq. 
(jA6|) . we obtain the motion equation 

OtPj(T ) = —i[H ay8 (T),pj(T)\ + £pj(T) (A12) 

for the generating density matrix pj(T ) with the initial 
condition pj(0) = p sys (0), where the Lindblad term is 

Cpj(T) = 2 E \/rv 70 i/ 9 j(r)oj cos[fc 0 (^ - z?)] 

ij 

- E UT”cos[fc 0 (xj - ^)]{OjOj,pj(T)}. (A13) 


Finally, by adding the Lindblad term describing the 
free space decay, we reproduce the results flUJ, (EH), and 
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